Sensitivity of neutron radii in a 208 Pb nucleus and a neutron star 
to nucleon-sigma-rho coupling corrections in relativistic mean 

field theory 

G. Shen, 1 J. Li, 1 G. C. Hillhouse, 1 - 2 <[] and J. Meng 1 - 3 ' 4 *! 

1 School of Physics, Peking University, Beijing 100871 
2 Department of Physics, University of Stellenbosch, Stellenbosch, South Africa 
3 Institute of Theoretical Physics, Chinese Academy of Science, Beijing 100080 
^Center of Theoretical Nuclear Physics, National Laboratory of 
Heavy Ion Accelerator, Lanzhou 730000 
(Dated: February 9, 2008) 

Abstract 

We study the sensitivity of the neutron skin thickness, S, in a 208 Pb nucleus to the addition of 
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relativistic mean field model. The PK1 and NL3 effective interactions lead to a minimum value of 
S = 0.16 fm in comparison with the original value of S = 0.28 fm. The S271 and Z271 effective 
interactions yield even smaller values of S = 0.11 fm, which are similar to those for nonrelativistic 
mean field models. A precise measurement of the neutron radius, and therefore S, in 208 Pb will 
place an important constraint on both relativistic and nonrelativistic mean field models. We also 
study the correlation between the radius of a 1.4 solar-mass neutron star and S. 
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I. INTRODUCTION 



Precise values for proton and neutron densities in nuclei, as well as their corresponding 
root-mean-square (rms) radii, are very important for providing quantitative predictions in 
nuclear physics and nuclear astrophysics. Our present understanding of nuclear phenomena 
largely hinges on accurate and model independent determinations of charge densities and 
radii via electron scattering. However, this picture of the nucleus will be more complete once 
precise values for neutron radii become available. Currently, the situation regarding the ex- 
perimental determination of neutron radii is unsatisfactory with errors typically being an 
order of magnitude larger than those for proton radii 1]. What is more disconcerting, how- 
ever, is the large variation in quoted values of the neutron radius for a single nucleus such as 

208p b 

0JMBQ 

. This dismal situation has prompted an experiment at Jefferson Lab- 
oratory [8| to measure the neutron radius in 208 Pb accurately and model independently via 
parity-violating electron scattering, the results of which are expected to have a widespread 
and important impact in nuclear physics. 

Reliable neutron densities and radii are needed for quantitatively characterizing the bulk 
properties of nuclear matter at normal nuclear densities, which in turn serve as valuable 
input for: studying the dynamics of heavy-ion collisions, determining the equation-of-state 
which underpins the evolution and structure of supernovae and neutron stars, studying the 
properties of exotic neutron-rich nuclei far from the beta-stability line, reducing uncertainties 



in atomic parity violation experiments,as well as 



or studying pionic and antiprotonic atoms. 



10, 



111 ] has been devoted to studying the 



Recently, much attention 
neutron skin thickness, 5*, in 208 Pb, a quantity which is defined as the difference between the 



rms radius of neutrons, ^/ < r 2 >, and protons, r 2 >, i.e., S = a/< r 2 > — < r 2 >. 
A precise measurement of the neutron radius of 208 Pb, together with existing high-precision 
measurements of the proton radius Q], will yield a precise value for the neutron skin thick- 
ness, thus providing one of the most stringent tests to date for current models of nuclear 
structure. Indeed, the Parity Radius Experiment at Jefferson Laboratory aims to measure 
the neutron radius in 208 Pb to an unprecedented accuracy of 1% (±0.05 fm). Currently 
relativistic and nonrelativistic models predict completely different values for the neutron ra- 
dius and, hence, it is evident that the latter measurement will have an impact on deepening 
our understanding of the dynamical basis of nuclear structure. In particular, nonrelativistic 
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Skyrme models predict S = 0.1 ~ 0.2 fm |9|, |l0j, |ll( , whereaspresent relativistic mean field 



(RMF) models, on the other hand, give S = 0.2 ~ 0.3 fm |9[ 11]. Model-dependent anal- 
yses of experimental data with hadronic probes yield values of S varying between 0.0 and 
0.2 fm (see Table EJ), which seem to be more consistent with the predictions of nonrelativistic 
models based on the Schrodinger equation. Hadronic measurements, however, suffer from 
potentially serious theoretical systematic errors associated with uncertainties in nuclear re- 
action mechanisms and, hence, one should be cautious about drawing conclusions regarding 
the appropriateness of various dynamical models. 

In this paper we focus on a relativistic description of nuclei and neutron stars within 
the field theoretical framework of quantum hydrodynamics. In particular, we capitalize 
on the success of the Walecka model and extensions thereof to successfully describe the 
properties of nuclear matter and finite nuclei 14, [3, 16, 17], as well as predict the 



well-known central and spin-orbit potentials which are usually postulated in nonrelativistic 
Schrodinger-equation-based models [3, 3, 0, Q • More specifically, we supplement exist- 
ing relativistic mean field Lagrangian densities with two new isospin dependent higher order 
correction terms relating to the coupling of the nucleon current to sigma- and rho-meson 
fields. These terms can be associated with multi-meson exchange processes occurring in the 
inner higher-density region of nuclei. In low-energy nucleon-nucleon scattering models, the 
contribution to the scattering amplitudes is usually dominated by single-meson exchange 
with an effective mass of less than 1 GeV. However, it has been demonstrated that the 
inclusion of two- meson exchange at low energies provides a major improvement in the de- 
scription of nucleon-nucleon scattering data 22]. In general there is no unique way for 
introducing nonlinear isovector-nucleon couplings in relativistic mean field models. In this 
paper we introduce the simplest structure generating the nonlinear isospin-dependent terms. 
Values for the various combinations of the new coupling constants are extracted by fitting 
to the properties of nuclear matter - such as saturation density, binding energy per nucleon, 
nuclear incompressibility and the symmetry energy - and the corresponding values for the 
neutron skin thickness in 208 Pb are extracted for a selection of RMF models. In particular, 
we focus on the PK1 



NL3 |24, S271 and Z271 [25( RMF parameterizations. 
One of the issues we wish to address in the paper is to find ways of adjusting the neutron 
radius in 208 Pb by adding new terms to the Lagrangian density of existing relativistic mean 
field models, while at the same time keeping the value of proton radius (which is accurately 
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determined experimentally) fixed. In particular, the question arises as to whether it is 
possible for RMF models to produce values of the neutron radius which are in the same 
range as those for nonrelativistic models. 

In recent analyses, a linear relationship between the neutron skin in 208 Pb and the nuclear 
symmetry energy at saturation density was proposed We also investigate the 



relationship between the latter quantities for our new RMF models. The new interaction 
terms considered in this paper influence the density dependence of symmetry energy, which 
in turn affects the neutron skin in 208 Pb. In addition, we extrapolate from normal to dense 



neutron matter and study the correlation between the radius o 
neutron radius in 208 Pb. Indeed, Horowitz and Piekarewicz 



a neutron star and the 
have performed such a 



correlation analysis for a wide range of RMF models and they concluded that, whereas the 
radius of a 0.5 solar-mass neutron star can be inferred from the radius of 208 Pb, the radius of 
a 1.4 solar-mass neutron star is not uniquely constrained by a measurement of the neutron 
skin thickness. In this paper we perform a similar analysis for the PK1, NL3, S271 and Z271 
effective interactions supplemented with our new terms, but for simplicity we only consider 
the mass and radius of a 1.4 solar- mass neutron star. 

This paper is organized in the following sequence. A brief formulation of our new RMF 
model in relation to existing RMF models for finite nuclei and neutron star matter is pre- 
sented in Sec. |H] The extraction of parameter sets for our new RMF models is discussed 
in Sec. IIHI In Sec. IIVI results are presented and discussed for the neutron radius in 208 Pb. 
Values of neutron radii in 208 Pb for different parameter sets are then correlated with the 
radius of a 1.4 solar-mass neutron star. Finally, in Sec. we summarize the main points 
of this paper. 

II. HIGHER ORDER CORRECTIONS TO RMF MODELS FOR SPHERICAL NU- 
CLEI AND NEUTRON STARS 



The basic physics und erly ing RMF models and their application in nuclear physics can 
be found in Refs. j3, Q, 0|- I n this section we present our new RMF model and indicate 
its relation to existing RMF models. The basic ansatz of the RMF theory is a Lagrangian 
density whereby nucleons are described as Dirac particles which interact via the exchange 
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of sigma- (a), omega- (u), and rho- (p) mesons, and also the photons (A), namely: 

,l + r 3 



C = V 



i^dft - m- g a a- gu'fu^ - g p ^f ■ p M - e7 M 
1„„ „ 1 



-A, 



- ~A» V A 
4 



(1) 



where the field tensors of the vector mesons and the electromagnetic field take the following 
form: 



d"u v - <9W 



A» v = d»A u - d u A p , 



(2) 



The nonlinear self-couplings for a and u mesons are respectively: 



U{a) = W + W, U{u>) = ±c 3 (u^ p f , 



(3) 



with the self-coupling constants g 2 , g 3 , and c 3 . 

To modi fy th e density dependence of nuclear matter symmetry energy, Horowitz and 
Piekarewicz [25] introduced the following nonlinear omega-rho coupling term: 



C H p = ^Kglpn ■ f^gluj^ 



(4) 



In addition to considering the above term, in this paper we introduce, for the first time, 
the following two new couplings of the nucleon current to sigma- and rho-meson fields: 



d = -Y 1 i>g p ^{g (T (T/m) f ■ /y0 



(5) 



C 2 = -T 2 ipg P Y(9acr/m) 2 T- p p tp 



(6) 



i.e., the total Lagrangian density of interest is: 



C'= C +C H p + £i + C 2 



(7) 
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The classical variational principle leads to the following equations of motion: 



[a-p + V(r) + p(m + 5(r))]-0i = e^i 



(8) 



(9) 



(10) 



for the nucleon spinors, where 

V(r) = (3{g^ + g p (l + V^a/m) + Y 2 (g a a /mf)r ■ fa + e^-±^ A^} 
5(r) = ^ CT (j 

and 

(ml - V 2 )(J= -gaps - g2° 2 - g 3 (J 3 - g P Pt,-J^(Ti(g a (r/m) + 2Y 2 (gla/m 2 )) 
(m 2 v - V 2 K= g u ji* - c 3 ^(u;^) - %Kg 2 p 9lh ' 

(m 2 - V 2 )/^ = 5,^(1 + iM&^/m) + r 2 (( ?(j( 7/m) 2 )-8A^ 2 ^^ 
- V 2 A^ = ejg 

for the mesons and photons. The nucleon spinors provide the relevant source terms: 

Ps = Etl Mi 

f = Ef=i ^il^i 
f = Ef=i ^il^i 

where the summations run over the valence nucleons only. The present method neglects 
the contribution of negative energy states, i.e. the so-called no sea approximation. The 
above non-linear equations are solved by iteration within the context of the mean field 
approximation whereby the meson field operators are replaced by their expectation values. 



A. Spherical nuclei 

Since spherical nuclei such as 208 Pb respect time reversal symmetry, there are no currents 
in the nucleus and the spatial vector components u>, p and A vanish. One is left with 
the time-like components uq, po and Aq. Charge conservation guarantees that only the 
3-components of the isovector p ,3 survive. As the system obeys rotational symmetry, the 
nucleon potentials and the meson-field sources depend only on the radial coordinate r. The 
Dirac spinor, ip(r) is characterized by the angular momentum quantum numbers l,j,m, 
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the isospin t = +(— ) for protons (neutrons) respectively, and any other relevant quantum 
numbers i, i.e. 

where Yj m (9,<f)) are the usual spinor spherical harmonics, G l /(r)/r and F-\r)/r denote 
the upper- and lower- component radial wave functions, and Xt represents the isospin wave 
functions specified by the isospin t. The radial wave functions are normalized as: 



/>oo 

/ dr(\G l /(r)\ 2 + |if(r)| 2 ) = l. (12) 
Jo 

The following phase convention is adopted for the vector spherical harmonics: 

(a-f)Y l = -Y 1 ' 

where 

l' = 2j-l = 



l + l, J = 1 + 1/2, 
I -I, j = l- 1/2. 

After some tedious algebra, one gets the radial equation for the upper and lower compo- 
nents, respectively jl^j ]: 



6,G?(r) = (-£ + f) if (r) + (m + S(r) + V(r))G l ^r) 

, = (+& + f)^'(r) - (m + S(r) - V(r))F?(r) 

where 

-(j + 1/2), for j = 1 + 1/2, 



(13) 



and 



(14) 



+U + 1/2), lor j = 1-1/2, 

V(r) = + g p (l + Y x {g a a/m) + T 2 (g e7 a/m) 2 )r 3 p 0!3 + \e(l + t 3 )A , 

S(r) = g a a. 

The meson fields equations reduce to radial Laplace equations of the form: 

~ r d~r + m V = ( } 
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where are the meson masses for <p = a, u, p and = for the photon. The relevant 
source terms are: 

-9aPs - 92<y 2 ~ 93V 3 ~ gpPo,3p3{Ti{ga(T/ni) + 1Y 2 {glajm 2 )) for the a field, 



S 6 = < 



g^Pv - c 3 ^q- ^Kglglpl^Q for the u field, 

(16) 



g p pz{\ + F^a/m) + ra^a/m) 2 ) - 8A,, g 2 p g^p 0t3 for the p field, 

ep c (r) for the Coulomb field, 

where the various nucleon densities are given by: 

4vrr 2 p s (r) = Eii(l^(r)| 2 - |F,(r)| 2 ), 

A^p v {r)= Ef=i(l^(r)| 2 + |F,(r)| 2 ), 

4vrr 2 p 3 (r) = Ef=i(l^(r)| 2 + |F,(r)| 2 ) - Ef=i(l^(r)| 2 + |F,(r)| 2 ), 

[4vrr 2 p c (r)= Ef =1 (l^(r)| 2 + |F,(r)| 2 ). 

The neutron and proton rms radii are directly related to the neutron and proton density 
distributions via the following relationship: 



where pj(r) (i = n,p) denotes the corresponding neutron and proton baryon density distri- 
butions. 

The total binding energy of the system is: 

E = E nudeon + E a + E p + E u + E c — mA 

= X/* - 2 / d3r {g^Ps(r) + -# 2 ^ 3 + 2^ 4 } 

i 

-\f d 3 rg pPo , 3 {(l + 2T 1 9 -^- + 3 r 2 (^) 2 )p 3 (r) - 8A,(^ ^po, 3 ) 2 } 
~\j d 3 r{ gu) uJo - ^c 3 u*}p v (r)- ^ J d 3 rA p c (r)-mA. (19) 

Eq. (13) is solved self-consistently using the Runge-Kutta algorithm and the shooting 
method, and the Laplace equations, given by Eq. (15), are solved employing a Green Function 
method. 

B. Neutron stars 

Next we discuss the RMF Lagrangian density used for describing a neutron star consisting 
of only nucleons and leptons (mainly electrons, e~, and muons, p~). The Lagrangian density 
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associated with the leptons is: 

^lepton 



^2 V'a(*7'*^/i - m \)^> 

X=e~ ,fi~ 



(20) 



such that the full Lagrangian density under consideration consists of a hadronic part, given 
by Eq. (J7|), plus a leptonic part, given by Eq. (|2U|) . Introducing the mean field and no- 
sea approximations, the equations of motion for baryons and mesons can be derived. The 
detailed formalism of which can be found in Ref. [32j. The Dirac equation represents the 
equation of motion associated with the leptons. Based on chemical equilibrium and charge 
neutral conditions, for a given total baryon density, the equations of motion are solved self- 
consistently by iteration, and the meson fields, particle densities and Fermi momenta of each 
species (protons, neutrons, electrons, and muons) are obtained simultaneously. 

The canonical energy-momentum tensor is derived by invoking the invariance of spacetime 
translation, namely: 

dC 



-g 



(21) 



where g^ u is the Minkowski metric in rectilinear coordinates. For the case in which neutron 
star matter can be considered to be a perfect fluid, the energy-momentum tensor is given 
by Q: 



<jpv 



(22) 



( e o o o\ 

p 
OOpO 
\0 p/ 

Comparing the ground-state expectation value of Eq. (|2*Tf with Eq. (J22J), via the field equa- 
tions the energy density, e, and the pressure, P, are given by: 
1 111 31 



+4E 

IT 2 ^ 

1 



7T In 

B=n,p J u 
1 

3 1 



k 2 dk\/k 2 + (mg 



r> ' 2 2 1 :i 1 [ 1 2 2 L i 1 2 2 i 2 2 \ 2 2 

P = -T> m a° ~ o9i° ~ 193° + -m,w + -c 3 uj + -m p p 03 + Ag p g w A v u oPQZ 



1 

-i 

2 



3tt 2 



E 

B=n,p 



fcfl 



o \/k 2 + [m B + g a o) 2 



dk 



7T 2 ^ 

A=e ,/i 

1 „ 1 

1 



k 2 dkJk 2 + m\, (23) 



3tt 2 



E 



\=e~ ,fi~ 



(24) 



where ks and k\ are the Fermi momenta of baryons and leptons respectively. 

The mass and radius of a neutron star are obtained by employing the Oppenheimer- 
Volkoff (OV) equations 

dp(r) [p(r) + e(r)] [M{r) + 4nr 3 p(r)] 
dr r(r - 2M(r)) ' ^ ' 

M(r) = 4tt f e(r')r' 2 dr', (26) 



where r denotes the radial coordinate relative to the center of the star, p(r) and e(r) are the 
pressure and energy density at a radial point r in the star respectively, and M(r) represents 
the mass of the sphere contained within a radius r. Since zero pressure cannot support a 
neutron star from collapsing, we define the radius, R, of the star as that radius at which 
the pressure is zero. The mass total mass, M(R), of the star is subsequently defined as the 
mass contained within a sphere of radius R. 

III. EXTRACTION OF PARAMETER SETS FOR NEW RMF MODELS 

In this section we present the procedure and criteria for extracting values for our new 
coupling constants Ti and 1^2 in Eqs. (JSJ) and (JBJ) respectively. In particular, we consider 
the addition of these new isopsin-dependent higher order correction terms to the PK1 123], 

n 

NL3 [24| . S271 and Z271 |25( RMF models. The parameters sets for these interactions are 



presented in Table ITT1 



The PK1 effective interaction 



We start with the recent PK1 effective interaction 
provides an excellent description of the properties of nuclear matter as well as for nuclei near 
and far from the valley of (3 stability. For this interaction, symmetric nuclear matter satu- 
rates at a Fermi momentum of 1.30 fm _1 with a binding energy per nucleon of —16.27 MeV 
and an incompressibility of K — 283 MeV. For comparison with the results of the PK1 
effective interaction, we also employ the NL3, S271, and Z271 effective interactions, which 
have been used in Refs. [27] to observe the effects of adding a nonlinear omega-rho cou- 
pling term, given by Eq. (J3J), on the neutron skin thickness, S, in 208 Pb and the radius of 
1.4 solar-mass neutron star. The NL3 effective interaction has also been used extensively 
to reproduce a variety of nuclear properties, such as binding energies, nuclear radii, nuclear 
density distribution, single particle spectra, etc. The NL3 and S271 effective interactions 
contain a sigma-meson self-coupling and Z271 includes both sigma- and omega-meson self- 
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couplings. These three interactions produce the following properties for symmetric nuclear 
matter: saturation at a Fermi momentum of 1.30 fm _1 , a binding energy per nucleon of 
—16.24 MeV, and an incompressibility of K — 271 MeV. 

Later in this section we will dicuss the procedure for calibrating values of the I\ and T 2 
couplings, in Eqs. (j3J) and (jHJ), with respect to the symmetry energy, a quantity which we 
now briefly discuss. The energy of asymmetric nuclear matter can be expanded around the 
energy of symmetric nuclear matter - where symmetric nuclear matter is characterized by 
identical proton and neutron densities, p p and p n respectively - as: 

fCM) = fOM = 0) + t 2 a sym (p) + 0(t 4 ), (27) 

where p = p n + p p = , k F is the the Fermi momentum, a sym (p) is the symmetry energy, 
and the asymmetry, t, is defined as: 

t = P^llL . (28) 

Pn + Pp 

The symmetry energy, a sym (p), describes how the energy of asymmetric nuclear matter 
changes with the asymmetry, t. ^From Eq. (27), and employing Eq. (19), one can extract 
the following expression for the symmetry energy of nuclear matter in RMF: 

a (a) - ( d ^ M ) - 4 + 9 * (29) 

b p 

where E*p = k 2 F + m* 2 , m* = m + g a a is the effective nucleon mass, and the effective 
rho-meson mass, m* is defined as: 



m 2 p + 8A^^> 2 o 
m p q _^ Yi^- + r 2 (— ) 2 ^ 2 ' 



f 2 



There are two terms in the symmetry energy. The first term in Eq. (I29j) represents the 
increase in the kinetic energy of the system when there is a relative displacement of the 
neutron and proton Fermi energies, and the second term in Eq. (129)1 represents the coupling 
between effective rho-meson and an isovector-vector current which no longer vanishes in 
asymmetric nuclear matter. 

We now describe our procedure for extracting values of g p , A v , and Ti, with T 2 = 0, 
in the Lagrangian density given by Eq. (j7jl. The symmetry energy at saturation density, 
corresponding to kf = 1.30 fm _1 , is not well constrained by the binding energy of nuclei. 
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However, some average of the symmetry energy at full density and the surface energy is 
constrained by the binding energy . We adopt the procedure of Horowitz and Piekarewicz, 
whereby we choose values of A„ and I\ with T 2 = 0, and then adjust the value g p such 
that the effective interaction has a fixed symmetry energy at an average density of kj = 
1.15 fm -1 (p = 0.10 fm -3 ) and the binding energy per nucleon in 208 Pb lies in the range 
\Ej A — (E/A) exp \ < 0.005 MeV, where (E/A) exp represents the experimental value. The 
symmetry energy at kf = 1.15 fm -1 is 26.08 MeV for the original PK1 effective interaction 
and 25.68 MeV for the original NL3, S271, and Z271 effective interactions, respectively. A 
similar procedure is followed for extracting parameter values for g p , A v , and T 2 , with Ti = 0. 



The experimental binding energy per nucleon, (E/A) exp , in 208 Pb is —7.868 MeV 35j]. As 
the correction for center of mass motion is neglected in the self-consistent RMF calculation, 
for comparative purposes an "equivalent" experimental (E/A) exp = —7.843 MeV for 208 Pb 
is used instead, which is estimated via the harmonic oscillator phenomenological formula 
E cm = 0.75 hwo with ftw = 41A" 1 / 3 Q. 

To avoid abnormal solutions at zero density for nuclear matter, A„ in Eq. (j3J) should not 
be negative. However, there is no such constraint on the values of Ti and T 2 in Eqs. (5) and 
(6) respectively. The role of the I\, T 2 and A„ couplings is to modify the density dependence 
of the symm etry energy. In Fig. Q we see that increasing Ti, decreasing T 2 , or increasing 
A„ (see Ref. [27]), causes the symmetry energy to grow more slowly with density. 

After readjusting the coupling constant g p , we perform self-consistent calculations as 
stated in Sec.|H]to determine the properties of 208 Pb and a 1.4 solar-mass neutron star. The 
numerical values of different combinations of A„, Ti (or T 2 ) and g p , and the corresponding 
binding energy per nucleon, E/A, the proton root mean square radius, R p , the neutron skin 
thickness, S, and the radius, R, of a 1.4 solar-mass neutron star, are also listed in Tables ITTT1 
to lVTTTl 

IV. RESULTS AND DISCUSSION 
A. Neutron radius of 208 Pb 

In this section we study the sensitivity of the neutron skin thickness, S, in 208 Pb with 
respect to the addition of two new terms, given by Eqs (5) and (6), to the original Lagrangian 
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densities for each of the PK1, NL3, S271 and Z271 effective interactions. In particular we 
also extract values of the neutron rms radius for RMF parameter sets which are consistent 
with experimental binding energies as well as with the properties of nuclear matter specified 
in the previous section. 

We start by considering the PK1 effective interaction and study the effect of different 
combinations of A„ and V±, with r 2 = 0, as well as different combinations of A„ and T\, with 
r 2 = - in the total Lagrangian density given by Eq. ((7j) - on the neutron skin thickness, 
neutron and proton rms radii, as well as neutron and proton density distributions. In order 
to understand the role of various coupling constants, we initially focus on the special case 
where T 2 = and consider the effect of different combinations of A„ and V\, which are 
chosen according to the procedure outlined in Sec. 11111 The value of g p is always chosen so 
as to reproduce the symmetry energy at an average density of kf = 1.15 fm -1 corresponding 
to a asm = 26.08 MeV. In Fig. the binding energy per nucleon, E/A, versus the neutron 
skin thickness, S, is presented. Different lines correspond to results with different Ti for 
fixed values of as indicated in the figure. The direction of the arrow next to the various 
lines corresponds to increasing values of F\. The open circles are the results for T\ = 0. The 
symbols retain their meanings in similar figures to follow. The effect of increasing Ti is to 
decrease S, while decreasing the absolute value of the binding energy. With the additional 
constraint that the experimental binding energy per nucleon, (E/A) exp , is —7.843 MeV 
(indicated by the dashed line in Fig. |2j), one can extract values of S which range from 0.159 
to 0.277 fm for different combinations of T\ and A v . More detailed results are listed in Table 

nm 

The accurately measured value of the proton rms radius, R p , of 208 Pb, namely 5.45 ± 
0.02 fm 0, 

places a stringent constraint on the choice of coupling constants. In Fig. El we 
also display the neutron and proton rms radii as functions of different combinations of A„ 
and Ti, with r 2 = 0. It is seen that the neutron rms radius is sensitive to different combi- 
nations of Ti and A„, whereas the proton rms radius is insensitive and in good agreement 
with the experimental value. In Fig. HI the neutron and proton density distributions have 
been computed for five different A„ and Ti combinations, which reproduce the experimental 
binding energy per nucleon for 208 Pb. While combinations of A^ and Ti with small values 
of S tolerate larger neutron central densities, the proton density distributions essentially 
remain unchanged in all cases. 
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Next, we consider the PK1 effective interaction and study the effect of different combi- 
nations of A„ and r 2 , for I\ = 0, on S. The binding energy per nucleon, E/A, versus S is 
shown in Fig. |5J The symbols retain the same meanings as in Fig. |2J except that the roles of 
T 2 and Ti are exchanged. Constraining the values of S via the experimental binding energy 
per nucleon, yields essentially the same range of values for S as was obtained in Fig. |21 The 
main difference between Figs. 121 and 03 lies in the fact that the values of r 2 are negative and 
the absolute values thereof are larger than Ti for fixed values of A„, and for the same values 
of the neutron skin thickness. More detailed results for different T 2 and combinations 
are listed in Table HVl Once again, we mention the fact that the proton rms radius is well 
constrained by the experimental data. 

Previous RMF models predicted values of S in 208 Pb which ranged from 0.2 to 0.3 fm 

h|. In our present analysis, we obtain values of S ranging from 0.159 to 0.277 fm, for 
different A„ and Fx (or T 2 ) combinations for the PK1 effective interaction, with the binding 
energy and proton rms radius both constrained by data. More specifically, our new models 
yield a minimum value of S = 0.159 fm for A„ = 0.035 and Ti = 0.2 (or T 2 = —0.55} in 
comparison with S = 0.16 ± 0.02 fm given by modern Skryme Hartree-Fock models [10 1. 

We now repeat the above analysis for the NL3 effective interaction. In Fig. |U1 the binding 
energy per nucleon, E/A, versus S is shown for different combinations of (a) A„ and Ti (left 
panel), with T 2 = 0, and (b) for different combinations of A„ and T 2 (right panel), with 
Y\ — 0. Constraining the values of S via the experimental binding energy per nucleon, 
yields a range of values for S similar to those associated with the PK1 effective interaction. 
The minimum value of S is 0.155 fm and this corresponds to the A^ = 0.035 and T\ = 0.25, 
with T 2 = 0, combination. More detailed data regarding the parameter sets and associated 
observables are listed in Tables IVl and IVTl 

Next, we study the effect of adding our new terms on two additional effective interactions, 
namely the S271 and Z271 parameterizations of RMF j^l^l. The original parameters for 
these interactions were extracted by fitting to the properties of nuclear matter and the 
experimentally determined proton radius of 208 Pb. The addition of the A„ term, given by 
Eq. Q), to the Lagrangian Eq. (^), reduces the values of the neutron skin thickness in 
208 Pb for both the S271 and Z271 effective interactions, while sacrificing good fits to the 
experimental binding energy per nucleon. The latter was the result of Ref. [27| . the points 
of which are indicated by open circles in Fig. The effect of different combinations [in the 
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Lagrangian given by Eq. (JJJ)] of A„ and Fx, with F2 = 0, and various combinations of A„ and 
T 2 , with Ti = 0, on S is displayed in the different panels in Fig. Constraining the values 
of S via the experimental binding energy per nucleon for the S271 effective interaction (see 
upper panels in Fig. |7J), yields values of S which range from 0.106 to 0.258 fm, for different 
combinations of Fx (r 2 ) and A„. A slight narrower range is obtained for the Z271 effective 
interaction, namely S = 0.134 ~ 0.241 fm: see lower panels in Fig. [HI The main difference 
between the S271 and Z271 effective interactions is that the signs of Ti and r 2 are opposite 
at points corresponding to the experimental binding energy per nucleon. More detailed data 
on S versus the binding energy per nucleon for different Ti and A^ combinations for S271 
and Z271 are presented in Tables IV 111 and IV1111 As is the case with the PK1 and NL3 
effective interactions, we have checked that the proton rms radius of 208 Pb agrees well with 
the experimental data. 

Previously it has been demonstrated that, for the original RMF models of interest to 
this paper, a linear relationship exists between the neutron skin thickness for 208 Pb and 
symmetry energy of nuclear matter at saturation density After the inclusion of 

our new terms, given by Eqs. © and (JBJ), we study this relation for each of PK1, NL3, S271 
and Z271 effective interactions: the results are presented in Fig. |HJ For this study we only 
take combinations of A^ and Ti (r 2 ) which reproduce both the experimental binding energy 
as well as the experimental proton rms radius of 208 Pb. Indeed, we also observe a linear 
relation between the neutron skin thickness, 5* (in fm), and the nuclear symmetry energy, 
a sym (in MeV), which can be expressed as: 

S « -0.283 + 0.015 a sym . (31) 

At the saturation point, the symmetry energy, a sym , for the various interactions, ranges from 
29 to 38 MeV and the neutron skin thickness, S, ranges from 0.14 to 0.28 fm. 

The addition of our new higher order correction terms, expressed via Eqs. (0) and ©, to 
existing RMF models can be associated with some kind of multi-meson exchange processes 
occurring in the inner (higher density) region of nuclei. In low-energy nucleon-nucleon 
scattering models, the exchange of mesons with an effective mass of less than 1 GeV seems 
to be the dominant contribution to the scattering amplitudes. After one-meson exchange, the 
inclusion of two-meson exchange at low energies has been demonstrated to provide a major 
improvement in the theoretical description of nucleon-nucleon scattering data 22]. In the 
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inner region of a heavy nucleus such as 208 Pb, one could imagine that multi-meson exchange 
processes, expressed via Eqs. (J3J) and (JHJ), could play an important role in modifying the short 
range behaviour of the nucleon-nucleon interaction, by softening the density dependence of 
the symmetry energy at high density (the inner part of nuclei) and thus leading to smaller 
values of neutron skin thickness S compared to previous interactions which excluded our 
new correction terms. 

In principle one could include the isospin-dependent nuclear force associated with the 
long-range pseudoscalar pion. However, because of parity conservation, the pion cannot 
contribute in a mean-field treatment of the ground state of spherical nuclei as 208 Pb. More- 
over, what determines the neutron skin thickness seems to be the neutron distribution in the 
inner core of nuclei, which in turn is more sensitive to the short-range part of the nuclear 
force. The above arguments suggest that pion contributions do not play an important role 
in determining the neutron skin thickness in 208 Pb. 

Since the parity- violating electron scattering experiment at Jefferson Lab can only provide 
information about the neutron radius, and hence also S, in 208 Pb, it is also important to 
study the consequences of an exact measurement of S on other properties of this nucleus. 
In particular, we now observe the influence of different Ti (T 2 ) and combinations (which 
reproduce both the experimental binding energy as well as the experimental proton rms 
radius of 208 Pb), for both PK1 and NL3 effective interactions, to neutron single particle 
energy levels versus their corresponding neutron rms radii in 208 Pb: the results are displayed 
in Fig. El We only consider the PK1 and NL3 effective interactions since they have been 
used to provide an excellent description of nuclei near and far from the valley of (3 stability. 
First we consider the original PK1 effective interaction (A„ = Ti = T 2 = 0) which yields S = 
0.277 fm. The corresponding neutron single particle energy levels are denoted by filled circles 
in the left panel of Fig. El We now compare these single particle energy levels to the values 
corresponding to those combinations of A v and T\, with T 2 = 0, and A v and T 2 , with Ti = 0, 
which yield a minimum value for S. In particular, the combination of A„ = 0.035 and Ti = 
0.2, with T 2 = 0, yields a mininum value of S = 0.159 fm for the PK1 effective interaction 
and the corresponding single particle energy levels are denoted by open squares in Fig. El 
Also the combination of A^ = 0.035 and T 2 = —0.55, with r x = 0, yields a mininum value of 
S = 0.159 fm: the corresponding single particle energy levels are denoted by open triangles 
in Fig. El The quantum numbers associated with each of the levels are the same as those 



16 



231 ] . The effect of decreasing 



associated with usual RMF predictions: see for example Ref. 
S is to deepen each single particle energy level and to decrease the corresponding neutron 
rms radius. More specifically, we see the addition of our new higher order correction terms 
is to deepen the deep-lying states more than the low-lying states. This, in turn, suggests 
the possibility that the effect of these new terms is to modify the short range behavior of 
the nucleon-nucleon interaction. In the right panel of Fig. El we observe the effect of adding 
our new terms on the neutron single particle energy levels and their corresponding neutron 
rms radii for the NL3 effective interaction: the results and discussion are similar to those 
associated with the PK1 effective interaction. 



B. Correlation between neutron skin thickness and the radius of a neutron star 

We now discuss the effect of our new terms, given by Eqs. © and (jOJ, on the radius, 
R, of a 1.4 solar-mass neutron star and also study its correlation with the neutron skin 
thickness, S, in 208 Pb for the PK1, NL3, S271 and Z271 effective interactions. In partciular, 
we consider the following combinations of parameters: (1) A„, I\ and T 2 terms separately, 
(2) A„ and T\ terms, with r 2 = 0, and(3) and r 2 terms, with T\ = 0. 

After adding the A„ term, given by Eq. to the Lagrangian given by Eq. the rela- 
tionship between R and S is displayed in Fig. EH different line types correspond to different 
interactions. For all the interactions, values of R increase with increasing values of S. How- 
ever, different interactions give different values of R for the same value of S. For example, 
for 5* = 0.21 fm, values of R vary from 13.27 km for the NL3 effective interaction to 11.63 km 
for the Z271 effective interaction. Therefore, the radius of a 1.4 solar-mass neutron star is 
not uniquely constrained by a measurement of the neutron skin thickness, as has already 
been pointed out by Horowitz and Piekarewicz [2jj . The reason for this is that S depends on 
the equation of state (EOS) near or below saturation density where most RMF models are 
calibrated to successfully describe the properties of nuclear matter and finite nuclei, whereas 
the neutron star radius is mainly sensitive to the EOS for higher densities. Hence, different 
RMF models give different predictions for the properties of neutron stars, and as such one 
cannot expect to reliably extrapolate any of these models to extreme conditions of isospin 
and density. However, combining separate measurements of S and R can provide valuable 
information on the EOS at low and high densities. For example, relatively large values of S 
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suggest a stiff low-density EOS, while relatively small values of R imply a soft high density 
EOS, and vice versa. 

In Fig. ^2 we observe the effect of the I\ term, given by Eq. (jSJ), with A^ = T 2 = 0, 
on the correlation between S and R. As expected, R increases with an increasing S, and 
different models give different relationships between R and S. The results for the PK1 and 
NL3 effective interactions are generally similar, but the value of R given by PK1 is smaller 
than that of NL3 for the same value of S. The S271 and Z271 effective interactions produce 
smaller neutron skin thicknesses, and also R increases more rapidly with S compared to 
PK1 and NL3. For example, for r x = 0, S ~ 0.280, 0.278, 0.254, 0.241 fm, and the 
corresponding values of R are 14.08, 13.89, 13.62, 13.07 km for the NL3, PK1, S271 and 
Z271 effective interactions, respectively. As the value of Ti is increased to 1.0, the neutron 
skin thickness decreases to 0.252, 0.253, 0.240, 0.236 fm, and the corresponding values of 
R decrease to 13.01, 12.74, 12.50, 11.97 km for the NL3, PK1, S271 and Z271 effective 
interactions, respectively. 

In Fig. El we observe the effect of the r 2 term, given by Eq. with A„ = Ti = 0, 
on the correlation between S and R: the results and discussion are very similar to those 
associated with Fig. [TT] 

The effect of different combinations of A„ with Ti, with r 2 = 0, and A„ with r 2 ,with 
Ti = 0, on the correlation between S and R for the (a) PK1, (b) NL3, (c) S271 and (d) Z271 
effective interactions is displayed in Fig. ED It is clear that the A^ term, given by Eq. (0J), 
plays the most prominent role for decreasing the values of R and S, i.e. the A„ term has 
the largest effect on changing values of R and S. In the upper left panel, associated with 
the PK1 effective interaction, the minimum neutron skin thickness in 208 Pb is 5* = 0.151 fm 
for the combination A„ = 0.02, Ti = 0.8, and T 2 = 0.0, while the minimum radius of a 1.4 
solar-mass neutron star is as low as R = 12.58 km when A^ = 0.01, Ti = 1.0, and T 2 = 0.0. 
The corresponding results for the NL3, S271, and Z271 effective interactions are shown in 
the other panels in Fig. ED The lowest value of the neutron skin thickness in 208 Pb is S = 
0.136 fm for the S271 effective interaction with A„ = 0.05, T\ = 0.0, and T 2 = —2.0, and 
the smallest value of R is 11.28 km for Z271 with A„ = 0.1, T 1 = 0.0, and T 2 = 0.0. 
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V. SUMMARY 



After introducing higher-order nucleon-sigma-rho coupling corrections to existing RMF 
models, and calibrating the various coupling constants with respect to the properties of 
nuclear matter, we observe a new range of values for the neutron skin thickness, S, in 208 Pb, 
namely 0.16 ~ 0.28 fm for the PK1 and NL3 effective interactions and 0.11 ~ 0.26 fm for 
the S271 and Z271 effective interactions. This means that, on one hand, for some particular 
effective interactions (S271 and Z271) it is possible to obtain values of the neutron radius in 
the same range as those for nonrelativistic mean field models, while on the other hand, for 
other effective interactions (PK1 and NL3) the neutron skin thickness is slightly extended 
compared to the previous upper limit of 0.2 fm. 

The addition of our new terms also has the effect of softening the density dependence 
of the symmetry energy at high densities. We also observe a linear relation between the 
neutron skin thickness and the nuclear symmetry energy at the saturation point. 

We have extrapolated our new relativistic mean field models from normal to dense neutron 
matter so as to correlate the radius of a 1.4 solar- mass neutron star with the neutron skin 
thickness in 208 Pb. We observe that the radius of a 1.4 solar- mass neutron star is not uniquely 
constrained by a measurement of the neutron skin thickness and is also model-dependent. 
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TABLE I: Summary of current values for the neutron skin thickness, S, in Pb 



Probe 


S (fm) 


Error (fm) 


Reference 


7r + and 7r~ scattering 


0.0 


0.1 


[2] 


Proton scattering (650 MeV) 


0.20 


0.04 


\3] 


Giant dipole resonance excitation 


0.19 


0.09 


\A] 


Nucleon scattering (40 - 200 MeV) 


0.17 




15} 


Proton-nucleus scattering (0.5-1.04 GeV) 


0.097 


0.014 


[6] 


Anti-protonic atoms 


0.15 


0.02 


m 



TABLE II: Parameter sets associated with the original PK1, NL3, S271, and Z271 RMF models. 
The nucleon and rho meson masses are kept fixed at M = 939 MeV for both protons, p, and 
neutrons, n, and m p = 763 MeV for the NL3, S271, and Z271 effective interactions. For the PK1 
effective interaction, M n = 939.5731 MeV, M p = 938.2796 MeV and m p = 763 MeV. 

RMF model g ff g u g 2 (fin" 1 ) g 3 c 3 m a (MeV) m w (MeV) 
PK1 10.3222 13.0131 -8.1688 -9.9976 55.636 514.0819 784.254 
NL3 10.217 12.868 -10.431 -28.885 508.194 782.5 
S271 9.006 10.806 -12.37 -17.323 505 783 

Z271 7.031 8.406 -5.435 63.691 49.941 465 783 
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TABLE III: Results for the PK1 effective interaction with T2 = 0. For each A„, the range of T\ is 
given with the corresponding g p coupling constant. Also listed are the binding energy per nucleon, 
E/A, the proton root mean square radius, R p , the neutron skin thickness, S, in 208 Pb, and the 
radius, R, of a 1.4 solar-mass neutron star. 



A„ 


Ti 4g 2 p E/A (MeV) R p (fm) S (fm) R (km) 


0.005 


0.00000 87.75942 -7.85382 5.44398 0.26332 13.5650 
0.05000 90.55426 -7.84285 5.44407 0.26224 13.4509 


0.01 


0.00000 93.54758 -7.86507 5.44561 0.24868 13.3450 
0.10000 100.40040 -7.84225 5.44608 0.24558 13.2510 


0.015 


0.00000 101.28410 -7.87067 5.44764 0.23500 13.1970 
0.15000 113.12450 -7.83869 5.44906 0.22817 13.0560 


0.02 


0.00000 109.87232 -7.87696 5.45031 0.22105 13.0880 
0.15000 123.96595 -7.84521 5.45239 0.21255 12.9596 


0.025 


0.00000 120.07776 -7.88201 5.45351 0.20705 13.0080 
0.20000 143.71214 -7.83903 5.45747 0.19295 12.8760 


0.03 


0.00000 132.38804 -7.88579 5.45726 0.19285 12.9440 
0.20000 161.69666 -7.84229 5.46243 0.17613 12.8290 


0.035 


0.00000 147.51000 -7.88824 5.46158 0.17828 12.8950 
0.20000 184.83000 -7.84363 5.46817 0.15861 12.7910 
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TABLE IV: Results for the PK1 effective interaction with Ti = 0. For each A„, the range of T2 is 
given with the corresponding g p coupling constant. Also listed are the binding energy per nucleon, 
E/A, the proton root mean square radius, R p , the neutron skin thickness, S, in 208 Pb, and the 
radius, R, of a 1.4 solar-mass neutron star. 



At, 


T 2 Ag 2 p E/A (MeV) R p (fm) S (fm) R (km) 


0.005 


0.00000 87.75942 -7.85382 5.44398 0.26332 13.5650 
-0.10000 89.34030 -7.84601 5.44358 0.26214 13.4870 


0.01 


0.00000 93.54758 -7.86507 5.44561 0.24868 13.3450 
-0.30000 99.48068 -7.84046 5.44480 0.24426 13.1470 


0.015 


0.00000 101.28410 -7.87067 5.44764 0.23500 13.1970 
-0.30000 107.86900 -7.84779 5.44734 0.22936 13.0240 


0.02 


0.00000 109.87232 -7.87696 5.45031 0.22105 13.0880 
-0.40000 120.51648 -7.84641 5.45060 0.21207 12.8950 


0.025 


0.00000 120.07776 -7.88201 5.45351 0.20705 13.0080 
-0.50000 136.42240 -7.84380 5.45487 0.19386 12.7980 


0.03 


0.00000 132.38804 -7.88579 5.45726 0.19285 12.9440 
-0.60000 157.20144 -7.83931 5.46024 0.17446 12.7290 


0.035 


0.00000 147.51000 -7.88824 5.46158 0.17828 12.8950 
-0.55000 175.83000 -7.84514 5.46548 0.15931 12.7150 
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TABLE V: Results for the NL3 effective interaction with 1^ = 0. For each A„, the range of Ti is 
given with the corresponding g p coupling constant. Also listed are the binding energy per nucleon, 
E/A, the proton root mean square radius, R p , the neutron skin thickness, S, in 208 Pb, and the 
radius, R, of a 1.4 solar-mass neutron star. 



A„ 


Ti 4g 2 p E/A (MeV) R p (fm) S (fm) R (km) 


0.005 


0.00000 84.86094 -7.86252 5.46039 0.26546 13.7800 
0.10000 90.36404 -7.84141 5.46047 0.26303 13.6790 


0.01 


0.00000 90.85902 -7.87088 5.46146 0.25128 13.5790 
0.15000 100.68116 -7.83939 5.46208 0.24585 13.4440 


0.015 


0.00000 97.85166 -7.87783 5.46320 0.23729 13.4390 
0.15000 109.28612 -7.84707 5.46441 0.23016 13.3270 


0.02 


0.00000 105.92526 -7.88395 5.46557 0.22333 13.3430 
0.20000 124.63490 -7.84268 5.46822 0.21123 13.2090 


0.025 


0.00000 115.51950 -7.88873 5.46853 0.20935 13.2670 
0.20000 138.10950 -7.84746 5.47226 0.19486 13.1530 


0.03 


0.00000 127.01290 -7.89232 5.47206 0.19521 13.2110 
0.25000 163.48180 -7.83926 0.17318 5.47856 13.0601 


0.035 


0.00000 141.00000 -7.89472 5.47619 0.18076 13.1630 
0.25000 187.42000 -7.84020 5.48457 0.15495 13.0500 
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TABLE VI: Results for the NL3 effective interaction with T\ = 0. For each A„, the range of I^ is 
given with the corresponding g p coupling constant. Also listed are the binding energy per nucleon, 
E/A, the proton root mean square radius, R p , the neutron skin thickness, S, in 208 Pb, and the 
radius, R, of a 1.4 solar-mass neutron star. 



At, 


T 2 Ag 2 p E/A (MeV) R p (fm) S (fm) R (km) 


0.005 


0.00000 84.86094 -7.86252 5.46039 0.26546 13.7800 
-0.30000 89.60516 -7.83969 5.45879 0.26154 13.5590 


0.01 


0.00000 90.85902 -7.87088 5.46146 0.25128 13.5790 
-0.40000 98.28740 -7.84082 5.45993 0.24456 13.3320 


0.015 


0.00000 97.85166 -7.87783 5.46320 0.23729 13.4390 
-0.50000 108.86836 -7.84069 5.46217 0.22699 13.1880 


0.02 


0.00000 105.92526 -7.88395 5.46557 0.22333 13.3430 
-0.60000 121.92576 -7.83948 5.46551 0.20865 13.0860 


0.025 


0.00000 115.51950 -7.88873 5.46853 0.20935 13.2670 
-0.60000 134.79210 -7.84435 5.46965 0.19256 13.0470 


0.03 


0.00000 127.01290 -7.89232 5.47206 0.19521 13.2110 
-0.70000 155.35130 -7.83993 5.47508 0.17258 12.9820 


0.035 


0.00000 141.00000 -7.89472 5.47619 0.18076 13.1630 
-0.65000 173.66000 -7.84539 5.48046 0.15725 12.9770 
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TABLE VII: Results for the S271 effective interaction with T2 = 0. For each A„, the range of Ti is 
given with the corresponding g p coupling constant. Also listed are the binding energy per nucleon, 
E/A, the proton root mean square radius, R p , the neutron skin thickness, S, in 208 Pb, and the 
radius, R, of a 1.4 solar-mass neutron star. 



A„ 


Ti 4g 2 p E/A (MeV) R p (fm) S (fm) R (km) 


0.01 


0.00000 91.58490 -7.94895 5.45986 0.23786 13.0880 
0.55000 121.352251 -7.84226 5.46037 0.22586 12.5867 


0.02 


0.00000 98.52548 -7.96073 5.46078 0.22151 12.7970 
0.65000 141.943385 -7.83858 5.46372 0.19862 12.3605 


0.03 


0.00000 106.66758 -7.97100 5.46247 0.20543 12.6290 
0.70000 165.27674 -7.84191 5.46869 0.17193 12.2640 


0.04 


0.00000 116.33780 -7.97978 5.46487 0.18950 12.5210 
0.75000 197.6836 -7.84262 5.47557 0.14352 12.3314 


0.05 


0.00000 127.87086 -7.98742 5.46795 0.17358 12.4470 
0.80000 245.799698 -7.83967 5.48465 0.11255 
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TABLE VIII: Results for the Z271 effective interaction with T2 = 0. For each A„, the range of I\ 
is given with corresponding coupling constant g p between rho meson and nucleon. Also listed are 
the binding energy per nucleon, E/A, the proton root mean square radius, R p , the neutron skin 
thickness, S, in 208 Pb, and the radius, R, of a 1.4 solar-mass neutron star. 



A„ 


Ti 4g 2 p E/A (MeV) R p (fm) S (fm) R (km) 


0.01 


0.00000 92.582884 -7.77718 5.45906 0.23478 12.6120 
-0.50000 79.602084 -7.83860 5.45937 0.23687 13.1190 


0.02 


0.00000 95.023504 -7.78256 5.45913 0.22827 12.2350 
-0.50000 81.396484 -7.84313 5.45928 0.23145 12.7610 


0.03 


0.00000 97.614400 -7.78766 5.45931 0.22182 11.9530 
-0.50000 83.283876 -7.84748 5.45928 0.22607 12.4490 


0.04 


0.00000 100.360324 -7.79255 5.45960 0.21542 11.7530 
-0.40000 88.021924 -7.84028 5.45937 0.21984 12.1030 


0.05 


0.00000 103.22560 -7.79745 5.46002 0.20901 11.6100 
-0.40000 90.212004 -7.84461 5.45960 0.21430 11.9140 


0.06 


0.00000 106.29610 -7.80200 5.46053 0.20267 11.5040 
-0.40000 92.544400 -7.84865 5.45990 0.20881 11.7690 


0.08 


0.00000 112.954384 -7.81084 5.46188 0.19004 11.3630 
-0.30000 101.123136 -7.84530 5.46099 0.19613 11.5130 


0.10 


0.00000 120.56040 -7.81876 5.46361 0.17751 11.2750 
-0.20000 111.386916 -7.84152 5.46268 0.18260 11.3530 
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FIG. 1: (Color online) The symmetry energy, a sym (in MeV), in symmetric nuclear matter as a 
function of the Fermi momentum, /cp (in fm -1 ), after adding different values of T\ (left panel) 
or T2 (right panel) terms to the original PK1 effective interaction. The direction of the arrows 
indicates values of increasing Ti and T2- 
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FIG. 2: (Color online) Binding energy per nucleon, E/A (in MeV), versus neutron skin thickness, 
S (in fm), in 208 Pb for the PKl effective interaction for different combinations of and T\, with 
r2 = 0. The direction of the arrow next to the various lines corresponds to increasing values of 
ri. The open circles denote the results for T\ = 0. The horizontal dashed line corresponds to the 
experimental binding energy per nucleon. 
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FIG. 3: (Color online) Proton and neutron radii, R p and R n respectively, in 208 Pb for the PK1 
effective interaction for different combinations of and Ti, with T2 = 0. The shaded regions 
denote the variation in the values of the radii associated with different combinations of these 
parameter sets. Also indicated, by the dashed horizontal lines, are the lower and upper limits 
associated with the experimentally extracted value of the proton radius. 
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FIG. 4: (Color online) Neutron, n, and proton, p, density distributions (in fm -3 ) in 208 Pb for the 
PK1 effective interaction for combinations of A,; and T\, with T2 = 0, which reproduce the the 
experimental binding energy per nucleon for 208 Pb. The value in square brackets indicates the 
neutron skin thickness, S, corresponding to a specific paramater set. 
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FIG. 5: (Color online) Binding energy per nucleon, E/A (in MeV), versus the neutron skin thick- 
ness, S (in fm), in 208 Pb for the PK1 effective interaction for different combinations of A„ and T2, 
with Ti = 0. The direction of the arrow next to the various lines corresponds to increasing values 
of FV The open circles denote the results for T2 = 0. The horizontal dashed line corresponds to 
the experimental binding energy per nucleon. 
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FIG. 6: (Color online) Binding energy per nucleon, E/A (in MeV), versus neutron skin thickness, 
S (in fm), in 208 Pb for the NL3 effective interaction for different combinations of (a) A v and F\ (left 
panel), with A2 = 0, and (b) and T2 (right panel), with A\ = 0. The direction of the arrows 
next to the various lines corresponds to increasing values of Ti (£2)- The open circles denote the 
results for Ti = (T2 = 0). The horizontal dashed line corresponds to the experimental binding 
energy per nucleon. 



34 



-7.76 


_ 1 1 1 1 1 1 1 J 


1 l\ 1 1 1 \ 1 - 


-7.8 


\ 0.02 
- 0.04 °-° 3 \ 
\0.05 \^ \ 


0.0l\ 0.00 \ - 


-7.84 






-7.88 






-7.92 






-7.96 

Q 


A v =0.06\ 

Nil I 1 


® (a) S271: A v + T, _ 





0.14 0.16 0.18 0.2 0.22 0.24 0.26 
S[fm] 



-7.8 

> -7.84 
S 

f -7.88 
= 

a -7.92 
-7.96 
-8 



_ I I I I I K 

\v0.03 

0.04 \. 
0.05 \ \ 


\ 0.02 


I i I i _ 
\ 0.01 \ 0.00 - 

\ \\ 






\\l 


- A v =0.06\ 






P 


I i 


\b) S271 : A, + T,_ 





0.096 0.116 0.136 0.156 0.176 0.196 0.216 0.236 0.256 
S[fm] 



-t 



-7.84 
-7.88 
-7.92 
-7.96 



I I I I I I I I I I 

o c 

o^ °°\ \ 


c 


i i i i _ 




i 


0.00 


r,\ \ \ \ \ o.oi 


\ \ \ 0.04 






\ \ 0.08 






0.12 






A,=0.15 (c)Z271:A v + r, 







III 



-7.8 - 



> -7.84 
I 

-7.88 



< 

a -7.92 



_ 1 1 1 1 1 1 1 1 

8 ° ° ° \ 
o^°°°\ \ 


1 


1 1 1 

1 : 


^\\ \ \ 

\ 0.08 
A v =0.16 0.12 




\ \ 1 0.00 
0.01 

04 

(d)Z271: A, + r 2 _ 





0.125 0.145 



0.165 0.185 0.205 0.225 0.245 0.265 0.285 
S[fm] 



0.125 0.145 0.165 0.185 0.205 0.225 0.245 0.265 0.285 
S[fm] 



FIG. 7: (Color online) Binding energy per nucleon, E/A (in MeV), versus neutron skin thickness, 
S (in fm), in 208 Pb for the (a) S271 effective interactions for different combinations of A„ and I\ 
(left panel), with A 2 = 0; (b) S271, A v and T 2 (right panel), with A 1 = 0; (c) Z271, A„ and Ti (left 
panel), with A 2 = 0; (d) Z271, and r 2 (right panel), with Ai = 0. The direction of the arrows 
next to the various lines corresponds to increasing values of Ti (r 2 ). The open circles denote the 
results for Ti = (r 2 = 0). The horizontal dashed line corresponds to the experimental binding 
energy per nucleon. 
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FIG. 8: (Color online) Neutron skin thickness, S (in fm), in 208 Pb versus the symmetry energy, 



''.syiii 



(in MeV), of nuclear matter at saturation density for the PK1, NL3, S271, and Z271 effective 



interactions for various combinations of and T± (T%), which reproduce the experimental binding 
energy per nucleon for 208 Pb. 
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FIG. 9: (Color online) Neutron single particle energy levels (in MeV) versus their corresponding 
neutron rms radii (in fm) in 208 Pb for the (a) PK1 and (b) NL3 effective interactions: original 
one (A,, = Ti = T2 = 0) denoted by filled circles, as well as for specific combinations of A„ and 
r\ (open squares), with T2 = 0, and also and T2 (open triangles), with Ti = which yield a 
minimum value for S (in square brackets) for the corresponding combinations of parameters. 



37 



a 



14.5 
14.0 - 
13.5 - 
13.0 - 
12.5 - 
12.0 
11.5 h 
11.0 




PKl 

NL3 
S271 
Z271 



0.13 0.16 0.19 0.22 0.25 0.28 

S[fm] 

FIG. 10: The radius, R (in km), of a 1.4 solar-mass neutron star versus the neutron skin thickness, 
S (in fm), in 208 Pb for the PKl, NL3, S271, and Z271 effective interactions for different ranges of 
values of A„, indicated next to the various line types. 



38 




0.23 0.24 0.25 0.26 0.27 0.28 0.29 

S[fm] 

FIG. 11: The radius, R (in km), of a 1.4 solar-mass neutron star versus the neutron skin thickness, 
S (in fm), in 208 Pb for the PK1, NL3, S271, and Z271 effective interactions for different ranges of 
values of T±, indicated next to the various line types. 
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FIG. 12: The radius, R (in km), of a 1.4 solar-mass neutron star versus the neutron skin thickness, 
S (in fm), in 208 Pb for the PK1, NL3, S271, and Z271 effective interactions for different ranges of 
values of T2, indicated next to the various line types. 
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FIG. 13: (Color online) The radius, R (in km), of a 1.4 solar-mass neutron star versus the neutron 
skin thickness, S (in fm), in 208 Pb for the (a) PK1, (b) NL3, (c) S271, and (d) Z271 effective 
interactions for various combinations of A„, T% and The upper panels are associated with the 
A„ and T\, with T2 = 0, combination: open circles denote values corresponding to T\ = and the 
direction of the arrows indicates values of increasing T\. The lower panels are associated with the 
A„ and T2, with Ti = 0, combination: open circles denote values corresponding to T2 = and the 
direction of the arrows indicates values of increasing IV 
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